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1.0  INTRODUCTION 


Free-space  range  specification  is  a  standard  method  used  to  assess  electromagnetic 
(EM)  system  performance.  The  specification  of  performance  of  a  free-space  range  parame¬ 
ter  enables  comparisons  of  different  systems  even  though  they  are  not  commonly  operated 
in  a  “free-space"  environment.  Free  space  propagation  may  be  defined  as  the  propagation 
of  energy  that  would  occur  if  an  omnidirectional  point-radiating  source  W'as  placed  in  outer 
space.  The  radiated  energy  would  travel  radially  outward  in  all  directions,  the  wave  fronts 
propagating  away  from  the  source  with  the  same  velocity  in  all  directions.  Obviously  these 
conditions  would  not  be  satisfied  if  the  point  source  was  placed  in  the  near-earth  environ¬ 
ment.  Refraction  by  the  atmosphere  insures  that  the  energy  is  not  propagated  with  the 
same  velocity  in  all  directions,  and  the  surface  of  the  earth  will  intercept  and  reflect  some 
portion  of  the  radiated  energy.  However,  the  concept  of  free-space  propagation  is  useful  for 
systems  operated  in  the  terrestrial  environment  because  such  propagation  can  be  used  as  a 
standard  measure  for  assessing  the  performance  of  these  systems.  The  free-space  detection 
range  of  a  radar  is  often  much  greater  than  the  actual  detection  range  of  that  radar  against 
targets  near  the  earth's  surface  because  of  horizon  effects.  Similarly,  the  free-space  intercept 
range  of  a  sensitive  Electronic  Support  Measure  (ESM)  receiver  may  be  several  hundred, 
or  several  thousand,  times  greater  than  the  actual  intercept  range  under  the  same  conditions. 

Utility  programs  have  been  developed  to  calculate  the  maximum  free-space  inter¬ 
cept  range  of  an  ESM  receiver  and  the  maximum  free-space  detection  range  of  a  pulse 
radar  system  against  a  specific-size  target.  These  utilities,  written  in  ANSI  Fortran,  are 
designed  to  be  used  as  simple  stand-alone  programs,  though  they  can  easily  be  incorpo¬ 
rated  into  other  programs  if  desired. 
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2.0  INPUTS,  OUTPUTS,  AND  LIMITS 


2.1  INPUTS 

The  determination  of  a  tree-space  maximum  range  requires  a  knowledge  of  certain 
EM  system  parameters.  The  number  of  parameters  required  depends  on  the  type  of  system, 
radar,  or  ESM.  More  parameters  are  required  for  a  radar  than  for  an  ESM  system.  The 
following  paragraphs  list  the  necessary  inputs,  their  units,  and  their  range  of  validity. 

2.1.1  Radar  System  Input  Parameters 

fable  2-1  lists  the  required  inputs  necessary  to  calculate  radar  free-space  detection 

range. 


Table  2-1.  Required  radar  and  target  parameters. 


Parameter 

I’nits 

Input  Range 

Simple  Radar 

n  a 

Yes  or  no 

Coherent  Integration 

n  a 

Yes  or  no 

Frequency 

MHz 

100.0  to  20.000.0 

Peak  Power 

kW 

>0.0  to  10.000.0 

Pulse  Width  (or  l.ength) 

MS 

>0.0  to  1.000.0 

Pulse  Repetition  Frequency 

Hz 

1.0  to  100,000.0 

Receiver  Noise  Figure 

dB 

>0.0  to  20.0 

Antenna  Gain 

dB 

>0.0  to  50.0 

Florizontal  Beam  Width 

deg 

>0.0  to  90.0 

Horizontal  Scan  Rate 

rpm 

1.0  to  500.0 

Hits  per  Scan 

integers 

1.0  to  1000.0 

Assumed  System  I.osses 

dB 

0.0  to  20.0 

Prooabilitv  of  Detection 

n  a 

0.1  to  0.9 

Prohabiltitv  of  False  Alarm 

n  a 

10  to  10  '  = 

farget  Radar  Cross-Section 

*> 

>0.0  lO  10’ 

Swerling  Case 

ii  a 

0  or  1 

Some  of  the  parameters  listed  may  not  be  required,  depending  on  the  type  of  radar  for 
w  hich  a  maximum  free-space  range  is  desired.  If  the  radar  is  a  simple,  scanning,  two- 
dimensional  (2D)  (i.e..  range  and  a/imuth)  system,  then  the  number  of  hits  per  scan  is  not 
entered  by  the  user.  In  this  case,  the  equivalent  variable,  the  number  of  pulses  integrated,  is 
calculated  f-'orn  the  horizontal  beam  width,  the  antenna  horizontal  scan  rate,  and  the  pulse- 
repetition-frequency  inputs.  Conversely,  for  more  complex  radar  systems,  such  as  height- 
finder  (.ID)  radars,  the  number  of  hits  per  scan  is  a  program  input.  In  this  case,  the  horiz¬ 
ontal  beam  width,  antenna  horizontal  scan  rate,  and  pulse  repetition  frequency  arc  not 
required,  and  the  number  of  hits  per  scan  must  be  entered  as  an  integer.  A  variation  of  the 
complex  radar  is  also  considered,  one  which  uses  coherent  integration,  the  best  example  of 
which  is  a  pulse-doppler  system,  fhe  user  is  required  to  specify  whether  the  radar  is  a 
simple  or  complex  system  and.  if  complex,  whether  or  not  coherent  integration  is  used  in 
processing.  The  entry  for  Swerling  Case  refers  to  the  type  of  target.  The  user  must  also 
specify  either  a  Swerling  Case  0  (steady,  nonductuating  target)  or  a  Swerling  Case  i  (slowly 
lluctuating  target ). 
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The  assumed-system-loss  parameter  is  used  to  include  all  miscellaneous  EM  system 
losses.  Such  losses  can  include,  but  are  not  limited  to.  collapsing  loss,  scanning  loss,  squint 
loss,  pulse-length  loss,  and  system-degradation  loss.  The  system-degradation  loss  occurs  as 
a  result  of  exposure  to  weather,  poor  maintenance  or  calibration,  and  deterioration  and 
aging  of  system  components.  A  value  of  5.0  dB.  or  more,  for  a  shipboard  radar  is  a  not 
uncommon  value  for  such  a  loss.  The  other  loss  mechanisms  are  described  by  Blake 
(1969,1980). 

2.1.2  ESM  Receiver  Input  Parameters 

Table  2-2  lists  the  EM  system  inputs  required  to  calculate  the  free-space  intercept 
range  of  an  ESM  receiver.  The  unit  abbreviation  dBm,  defined  as  decibels  above  1  mW, 
is  a  standard  way  to  specify  receiver  sensitivity  threshold  levels  for  ESM  systems.  The 
receiver  sensitivity  includes  the  gain  cf  the  receiving  antenna. 

Table  2-2.  Req  uired  E.SM  system  inputs. 


Parameter 

Units 

Input  Range 

Frequency 

MH/ 

100.0  to  20.000.0 

Transmitter  Power 

kW 

0.1  to  10.000.0 

Transmitter  Antenna  Gain 

dB 

5.0  to  50.0 

Receiver  Sensitivity 

dBm 

10.0  to  -150.0 

2.2  OUTPUTS 

The  only  outputs  are  the  calculated  radar  or  ESM  system  free-space  range  values 
in  kilometers  and  the  corresponding  one-way  free-space  path  loss  at  that  free-space  range. 
Sample  program  outputs,  and  their  corresponding  inputs,  will  be  given  in  section  4.0. 

2.3  LIMITS 

The  calculation  of  free-space  range  described  in  this  document  is  limited  to  pulse 
radars  for  operational  parameters  within  the  range  of  validity  of  the  inputs  of  Table  2-1 
and  ESM  systems  with  operational  parameters  within  the  range  of  validity  of  the  inputs  of 
Table  2-2. 


3 


3.0  MODELS 


There  are  two  different  models  used  to  determine  free-space  EM  system  range, 
depending  on  whether  or  not  the  system  is  a  radar  or  an  ESM  system.  The  model  for  the 
radar  system  will  be  discussed  first,  followed  by  a  discussion  of  the  model  for  the  ESM 
system. 

3.1  RADAR  FREE-SPACE  RANGE  MODEL 

Radar  free-space  range  determination  is  probabilistic,  since  the  detection  of  a  target 
depends  on  a  great  many  factors.  First,  the  target  must  be  in  the  radar  beam,  and  it  must 
be  there  long  enough  and  with  a  strong  enough  return  signal  that  it  can  be  recognized  as  a 
target  by  the  detector,  whether  or  not  the  detector  is  automatic  or  human.  Second,  the 
detector  must  be  able  to  discriminate  between  a  target  and  extraneous  background  noise.  If 
there  are  a  number  of  objects  at  the  same  range  as  the  target,  the  return  signal  from  these 
objects  may  mask  the  signal  returned  from  the  target.  These  objects  can  be  natural  sources, 
such  as  terrain  clutter  or  sea  clutter,  or  man-made  sources,  such  as  jammers  or  chaff.  Target 
return  signal  is  also  usually  a  complex  function  of  aspect  angle  and  may  vary  rapidly  as  a 
function  of  target  motion.  Atmospheric  refraction,  atmospheric  noise  from  extraterrestrial 
sources,  and  reflections  from  the  earth's  surface,  especially  the  ocean,  are  examples  of  other 
possible  natural  contributors  that  can  affect  the  detection  process.  Given  this,  it  should  be 
apparent  that  a  prediction  of  maximum  radar  range  cannot  be  guaranteed  to  be  accurate  in 
a  strict  sense.  Radar  performance  is  often  specified  in  terms  of  the  free-space  detection 
range,  which  simply  ignores  the  complicating  factors  which  arise  from  the  usual  opera¬ 
tional  environment  of  the  radar  system  and  concentrates  on  those  parts  of  the  detection 
process  that  are  free  of  this  contamination.  This  process  is  still  probabilistic,  because  any 
return  signal  will  still  be  mixed  in  with  noise,  but  it  provides  a  useful  reference  point  to 
radar  system  users. 

The  basic  maximum  radar  range  equation,  derived  from  Blake  (1980,  Eq.  1.28),  is 
given  by 


=  58.0 


Gr  p,  O  T 


fMH:'-  .V,  (5/.V),,„  L 


(1) 


w  here 


R,. 

P, 


O 


fMHz  ^ 


L 


maximum  radar  free-space  detection  range  in  kilometers 

transmitted  power  in  kilowatts 

radar  antenna  power  gain  ratio 

radar  target  cross-section  in  square  meters 

pulse  width  (or  length)  in  microseconds 

radar  system  frequency  in  megahertz 

minimum  signal-to-noise  ratio  for  a  specified  probability  of 
detection,  probability  of  false  alarms,  and  Swerling  Case  0  or  I 

receiving  system  noise  figure 

assumed  system  losses  expressed  as  a  ratio 
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This  equation  assumes  that  the  system  is  a  monostatic  radar  and  that  the  radar  receiver 
is  bandwidth-matched  to  the  pulse  width.  Equation  1  is  applicable  to  pulse  compression 
radars  when  r  is  the  width  of  the  compressed  pulse. 

As  previously  noted,  the  detection  process  is  statistical  or  probabilistic,  because  the 
returned  energy  from  the  target  is  mixed  in  with  receiver  noise.  Since  this  noise  is  random, 
there  will  always  be  voltage  fluctuations  at  the  output  of  the  detector.  These  fluctuations 
can  make  it  difficult  to  determine  with  any  certainty  whether  or  not  a  target  is  present. 
Rapid  voltage  changes  can  be  due  to  noise  or  target  return.  The  probability  that  a  target¬ 
generated  signal,  when  present,  will  be  detected  is  called  the  probability  of  detection.  Pj. 
and  the  probability  that  a  noise  fluctuation  will  be  mistaken  for  a  target  is  called  the  prob¬ 
ability  of  false  alarm,  The  minimum  required  signal-to-noisc  ratio  can  be  defined  in 

terms  of  these  quantities  by  using  an  empirically  derived  formula  (Blake.  1980.  Eq.  2.29); 

=  [■V„/(4  .V,,)]  M  Ml  +  <16  .V,,/.VJ]'  l  L,  (2) 

where 


A'„  =  (f(la  ^ 

(3) 

-  2.36  (  Log  Pfj'- 

1.02 

(4) 

II 

1^ 

(5) 

/  =0.9(2P,/-1) 

(6) 

Equation  2  is  sometimes  referred  to  as  the  detectability  factor.  .V^,  is  the  number  of  pulses 
integrated  by  the  detector  for  a  simple  2D  scanning  system.  If  the  radar  is  a  complex 
system,  height-finder  (.ID)  radar.  is  the  number  of  hits  per  scan.  .A'^,  for  simple  radars  is 
given  by  Blake  (1969.  Eq.  30): 

A’/,  =  «■>/, /Ml  A,,  >1.0  (7) 

where  (■);,  is  the  horizontal  beam  width  in  degrees. /rr/" is  the  pulse  repetition  frequency  in 
Hertz  and  hsr  is  the  horizontal  scan  rate  in  rpm.  Unlike  the  user-entered  number  of  hits  per 
scan.  .V^,  can  have  noninteger  values  but  is  limited  to  a  minimum  of  a  single  pulse  inte¬ 
grated.  Lj  is  the  lluct nation  loss  ratio  (/,,  =  1.0  for  a  steady.  Swerling  Case  0  target).  If  the 
target  is  fluctuating  (Swerling  Case  I ).  the  additional  loss  due  to  this  fluctuation  must  be 
calculated.  The  fluctuation  loss  ratio.  1.,.  is  given  by  Blake  ( 1980.  Eq.  2.45): 

L,  -  [(  In  PJ  ( 1  +  g^„)]  '  (8) 

Equation  2  is  valid  over  the  range  0. 1  <  /*,/  <  0.9  and  10  <  Pfa—  *0  '*■  this 

formulation,  a  square-law  detector,  a  uniform-weight  integrator,  and  a  constant  signal 
power  are  assumed.  Blake  ( 1980.  pp.  45  48)  notes  that  square-law  detectors  are  rarely  used 
in  radars,  but  the  difference  between  square-law  detectors  and  linear  detectors,  which  are 
more  commonly  used,  is  gcrierallv  much  less  than  1  dB.  While  Eq.  2  is  generally  valid  for 
simpler  radars  that  do  not  use  coherent  integration.  Blake  ( 1980.  p.  59)  notes  that  the  best 
wav  to  process  nonfluctuating  or  slowly  fluctuating  multiple-pulse  signals  is  to  integrate 
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them  coherently.  Such  processing,  he  notes,  is  achieved  only  in  pulse-doppler  radar  systems, 
though  some  other  forms  of  coherent  processing  are  used,  such  as  moving-target-indicator 
and  pulse-compression  svstems.  \\’hen  coherent  integration  is  used.  Eq.  2  becomes 


(•V  V),,,,,.  =  f  V„  (4  .V,,)lll  .[1  +(K'  .V,,)]'  !  /  ,  (9) 

where  all  quantities  are  as  previinisK  delined.  Equation  9  reveals  that  coherent  integration 
ot  .N'^,  pulses  is  equivalent  to  the  detectabilitv  factor  of  a  single  pulse  .V^,  times  as  long. 

Quantities  like  the  noise  figure.  .V,,  antenna  gam.  (i.  and  assumed  svstem  losses.  /.. 
ol  a  radar  system  are  usuallv  stated  in  decibels.  Since  f'q.  I  is  written  with  these  quantities 
m  terms  of  ratios,  it  is  convenient  to  restate  I  q.  1  with  these  terms  in  decibels.  Eo  accomp¬ 
lish  this,  another  quantity  is  defined  which  adds  all  the  previously  mentioned  terms  in 
decibels  and  allows  Eq.  1  to  be  given  in  terms  of  this  new  quantity; 

=  10.0  . . .  (10) 


Here  the  detectabilitv  lactor  Irom  liq.  or  I  q.  9.  whichever  applies,  is  also 
includetl  in  decibels  [  10  Eog(.S  M,,,,,,!  II  this  representation  is  used.  Eq.  1  becomes 


=  .‘iX.O 


fM  Hz- 


(km) 


(M) 


with  all  other  variables  retaining  their  previously  defined  units.  Eq.  1  1  is  the  most  conve¬ 
nient  e.xpression  lor  the  calculation  of  ladar  free-space  range  and  is  the  form  actually  used 
m  the  f  ortran  code 

Ev|uations  2  through  I  1  are  incorporated  into  a  Eortran  program.  RESDR  (Radar 
f  ree-Space  Detection  Range),  consisting  of  a  main  routine  and  three  subroutines.  When 
R  f  SDR  is  c.xecutcd,  the  main  routine  calls  two  input  subroutines.  INFI.XDR  and  IN  I  R(i  1 . 
to  prompt  the  operator  to  enter  the  required  radar  and  target  parameters  ol  lable  2-1. 

I'pon  return  Irom  the  INR.ADR  and  IN  I  R(i  I  subroutines,  the  free-space  range  is  calcu¬ 
lated  from  I  (.(  Il  l  he  equivalent  one-way  p;ith  loss  at  the  free-space  range,  assuming 
isotropic  antennas,  is  also  calculated  ;it  this  time  by  using  the  basic  transmission  equation 
from  Kerr(l9.Sl.  [  q  2-1,*'): 

I  OSS  =  .^2.4.S  +  200  I  o^(  f  Mf/z)  (12) 

The  free-space  loss  value  is  uselul  for  evaluating  radar  system  performance  with  trans- 
missiim-loss  plots  tor  tictual  operational  environments.  I  he  third  subroutine.  RDRDl'  E. 
is  then  called  to  tlisplav  a  list  of  the  input  parameters  and  the  calculated  free-space  range. 
Program  opertition  is  terminated  after  the  radar  and  tjirget  inputs,  the  calculated  free-space 
range,  and  the  path-loss  values  are  listed.  Complete  program  listings  for  RESDR  are 
included  in  Appendix  A. 

lo  use  RESDR.  the  operator  must  compile  and  link  the  F  ortran  routines  listed  in 
Appendix  A.  Numerictil  radar  system  or  target  input'  outside  the  specified  input  range  will 
be  program-limited  to  the  nearest  bound  specified  m  lable  2-1.  I  hat  is.  it  a  frequenev  of 
.‘vO  Mil/  IS  entered,  the  program  will  automaticallv  set  the  irequenev  variable  to  100  MH/. 
the  lower  frequenev  bound.  Blank  field  entries  are  not  acceptable  numerical  inputs,  i  hose 
inputs  listed  in  lable  2-1  thiit  require  ;i  numerictil  value  greater  than  zero  will  be  defaulted 
to  0  I  if  the  user  inputs  a  negative  or  zero  value.  I  his  is  to  prevent  runtime  errors  at 
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program  execution.  Non-numeric  (alpha-character)  inputs,  such  as  yes  no  responses  (e.g., 
type  of  system),  will  accept  any  character,  even  blanks.  The  default  for  such  an  entry  is  the 
first  option  listed  in  the  operator  prompt.  An  example  would  be  the  choice  of  the  type  of 
radar  system.  The  INTRGl  operator  prompt  asks:  “Is  the  radar  a  simple  system  [range 
and  a/imuth  only,  nonheight  finder,  simple  signal  processing]  (y  or  n)?”  Any  response, 
other  than  an  "n.”  will  cause  the  system  to  consider  the  input  as  a  simple  radar  system. 
Program  operation  is  as  stated  in  the  previous  paragraph. 

3.2  ESM  RECEIVER  FREE-SPACE  RANGE  MODEL 

I  hc  free-space  range  e'Hiation  for  an  ESM  receiver,  derived  from  Kerr  { 1951, 

Eq.  2- 1 5). is  given  by 

Rf  -  I  0  ()l  I  /’o  20  I  iiuO  U//rl  ♦  Of  S  •  27  S5I7)  20  0  ^  (13) 

where 

Ri^  =  free-space  intercept  range  in  kilometers 

P,  -  the  transmitted  power  in  kilowatts 

G,  -  gain  in  dB  above  an  isotropic  radiator  of  the  transmitter  antenna 

/'.V///r  =  frequency  in  megahertz 

S  -  sensitivity  of  the  receiver,  including  receiving  antenna  gain,  in  dBm 

The  range  values  obtained  by  using  Eq.  13  can  be  extremely  large,  perhaps  tens  of 
thousands  of  kilometers.  The  free-space  range  can  be  stated  as  an  equivalent  free-spt-ce 
loss  value  for  isotropic  antennas  by  using  Eq.  12.  This  one-way  path-loss  value  can  be  used 
to  evaluate  actual  ESM  receiver  intercept-range  performance  when  used  with  transmission- 
loss  plots  or  EM  propagation  codes.  The  ESM  free-space  range  program  returns  both 
the  range  from  Eq.  13  and  the  equivalent  free-space  path  loss  for  isotropic  antennas  from 
Eq.  12. 

Equations  13  and  12  are  incorporated  into  a  Fortran  program.  FSIR  (Free-Space 
Intei  cept  Range),  consisting  of  a  main  routine  and  two  subroutines.  When  FSIR  is  executed, 
the  main  routine  calls  INPUTS  to  prompt  the  operator  to  enter  the  necessary  EM  system 
parameters  of  1  able  2-2.  Upon  return  from  INPUTS,  ths  ;nain  routine  calculates  the  free- 
space  range  and  the  one-way  path  loss  at  the  free-space  range  by  using  Eq.  13  and  12, 
respectively.  The  main  routine  then  calls  ESMOUT.  This  subroutine  is  used  to  display  the 
input  parameters  and  the  calculated  free-space  range  and  path-loss  values.  Program  execu¬ 
tion  is  then  terminated.  Complete  program  listings  for  FSIR  are  included  in  Appendix  B. 

To  use  FSIR.  the  operator  must  compile  and  link  the  Fortran  routines  listed  in 
•Appendix  B.  Numerical  EM  system  inputs  outside  the  specified  input  range  will  be  program- 
limited  to  the  nearest  bound  specified  in  Table  2-2.  That  is,  if  a  frequency  ui  21.000  MEIz  is 
entered,  the  program  will  automatically  set  the  frequency  variable  to  20.000  MHz,  the  upper 
Irequency  bound.  Blank  field  entries  are  not  acceptable  numerica'  inputs.  The  program 
operational  sequence  is  detailed  in  the  preceding  paragraph. 
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4.0  TEST  CASES 


This  section  describes  test  and  evaluation  criteria  for  the  RFSDR  and  FSIR  pro¬ 
grams.  The  test  cases  will  allow  program  users  to  verify  correct  operation. 

4.1  RFSDR  TEST  CASES 

Table  4-1  lists  five  test  cases  that  should  verify  that  the  RFSDR  program  is  oper¬ 
ating  correctly.  The  input  and  output  parameters  are  listed  to  the  nearest  tenth  of  a  unit. 
RFSDR  is  operating  correctly  if  the  free-space  range  is  within  0.1  km  of  the  value  listed  in 
Table  4-1.  The  abbreviations  for  the  variables  are  the  same  as  those  listed  in  section  3. 1. 
with  the  following  exceptions.  The  radar  type  is  listed  as  either  an  “S"  (simple  2D  system) 
or  a  “C”  (complex  system).  If  the  radar  is  not  a  simple  system  it  may  use  coherent  process¬ 
ing.  In  this  case  the  line  labeled  “Coherent"  will  have  a  yes  or  no  entry.  The  target  Swerling 
Case  is  given  in  the  heading  labeled  “SW  Case."  The  one-way  free-space  path  loss  from 
Eq.  15  is  given  in  the  heading  labeled  “Toss."  An  entry  of  “n  a"  indicates  that  this  variable 
is  not  used  in  the  test  case. 

Table  4-1.  RFSDR  test  case  inputs  and  outputs. 


—  1 

Parameter 

Test  Case 

1 

3 

4 

5 

Radar 

S 

s 

S 

C 

C 

Coherent 

n  a 

n  a 

n  a 

No 

Yes 

tMH: 

100.0 

20.000.0 

100.0 

5000.0 

10.000.0 

P,  (A  H') 

1.0 

10.000.0 

10.000.0 

5000.0 

4000.0 

r  (AiV) 

0.1 

1000.0 

1000.0 

100.0 

1000.0 

N;  (dB) 

0  1 

20.0 

20.0 

20.0 

20,0 

(/  (dB) 

0.1 

50.0 

50.0 

.to.o 

35.0 

/.  (dB) 

0.0 

0.0 

0.0 

10.0 

5.0 

PrfiH/] 

1.0 

100.000.0 

100.000.0 

n  a 

n  a 

(■)^,  (deg) 

1.0 

90.0 

20.0 

n  a 

n  a 

Hsr  (rpm) 

1.0 

500.0 

500.0 

n  a 

n  a 

H'ts  Sean 

n  a 

n  a 

n  a 

80.0 

100.0 

a  (  m  - 1 

1.0 

1.0 

1.0 

10’ 

1.0 

Pj 

0.9 

0.9 

0.9 

0.1 

0.9 

10  I- 

10  •• 

10  1- 

10 

10 

SW  C  ase 

1 

0 

0 

0 

1 

( km ) 

4176.8 

42274.0 

7353,7 

346. 1 

1  OSS  (dB) 

190.9 

165,0 

183.8 

163.2 

4.2  FSIR  TEST  c  As  ' 

Table  4-2  lists  five  test  cases  that  should  verify  that  the  FSIR  program  is  operating 
correctly.  I  hc  input  and  output  parameters  are  listed  to  the  nearest  tenth  of  a  unit.  FSIR  is 
operating  correctly  if  the  free-space  range  is  within  0. 1  km  of  the  value  listed  in  Table  4-2 
for  test  cases  I.  2.  3.  and  5.  For  test  case  4.  the  number  of  significant  digits  precludes  using 
the  free-space  range.  In  this  test,  the  one-way  path  loss  should  be  within  0.1  dB.  The  abbre¬ 
viations  lor  the  variables  are  the  same  as  those  listed  in  section  3.2. 
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Tuhip  4-2.  FSIR  test  case  inputs  and  outputs. 


Parameter 

Test  Case 

1 

2 

3 

4 

5 

fMH: 

100,0 

100.0 

20.000.0 

20.000.0 

5500.0 

P,  (kW) 

0.1 

10,000.0 

10.000.0 

10,000.0 

1000.0 

G,  (dB) 

5.0 

50  0 

5.0 

50.0 

35.0 

S  (dBm) 

10,0 

10.0 

150.0 

150.0 

90.0 

Rf^  (km) 

0.4 

2.T855..1 

2.121.068.8 

1.193e9 

243.905.9 

I.oss  (dB) 

65.0 

160.0 

245.0 

300.0 

215.0 
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Appendix  A 

RADAR  FREE-SPACE  RANGE  (RFSDR)  PROGRAM  LISTS 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 


c 

c 

c 

c 


RFSDR  calculates  the  free-space  range  of  a  radar  and  the  pathloss 
at  the  free-space  range. 


VARIABLES : 

antgn 

dbterm 

floss 

fitihz 

hitscn 


hbwidth 

hscnrS 

icoher 

icmplx 

itgt 

pkpwr 

plsrS 

psubd 

psubfa 

pthlos 

rcvrnf 

rsubfs 

sigma 

snmin 

syslos 

tau 


DESCRIPTION: 

Radar  antenna  gain  in  dB 
Various  system  gains  and  losses  in  dB 
Fluctuating  target  loss  for  Swerling  Case  1  Target 
Radar  frequency  in  MHz 

Hits  per  scan,  specified  for  3D  radars  and 
calculated  for  2D  radars  as  the  number  of 
pulses  integrated 
Horizontal  beam  width  in  degrees 
Horizontal  scan  rate  in  rpm 

Coherent  integration  flag  [1  =  yes,  0  =  no] 

Complex  system  (3D  radar)  flag  [1  =  yes,  0  =  no] 
Fluctuating  target  flag  [1  =  fluct.,  0  =  steady] 
Peak  radar  power  in  kw 

Pulse  repetition  frequency  (or  rate)  in  Hz 

Probability  of  detection 

Probability  of  false  alarms 

Pathloss  at  radar  free-space  range  in  dB 

Radar  receiver  noise  figure  in  dB 

Radar  free-space  range  in  km 

Radar  cross-section  of  target  in  square  meters 
Minimum  signal-to-noise  ratio  in  dB 
Assumed  radar  system  losses  in  dB 
Pulse  width  (or  length)  in  microseconds 


Call  subroutine  incalc  to  enter  radar  parameters. 

real *4  antgn, dbterm, floss , fmhz ,gsubd,gsubf a, hitscn, hbwidth, hscnrS 
real*4  pkpwr , plsrS , pthlos , psubd , pusbf a , rcvrnf , rsubfs , sigma , snmin 
real*4  syslos, t, tau, XsubO 
inLeger*2  icoher , icmplx, itgt 


Call  subroutine  inradr  to  enter  radar  parameters. 

call  inradr (antgn, fmhz , hitscn, hbwidth, hscnrS , icoher, icmplx, 
1  pkpwr, plsrS , rcvrnf , syslos, tau) 

Call  subroutine  intrgt  to  enter  target  parameters. 

call  intrgt ( itgt , psubd , psubfa , sigma) 

T  =  .9*(2.*psubd  -  1.) 

gsubd  =  1.231*T  /  SQRT ( 1 . 0  -  T**2) 

gsubfa  =  2.36  *  SQRT (-LOG10 (psubfa) )  -  1.02 

XsubO  =  (gsubfa  +  gsubd) **2 

IF  ((icoher  .EQ.  1)  .AND.  (icmplx  .EQ.  1))  THEN 
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c  If  coherent  processing  is  tised  [pulse  doppler  radar]  icoher=l 

snmin  =  XsubO/4 . 0/hitscn* ( 1 . 0  +  SQRT(1.0  +  16.0/Xsub0)) 

ELSE 

c  Incoherent  integration. 

snmin  =  XsubO/4 , 0/hitscn* ( 1 . 0  +  SQRT(1.0  +  16 . 0*hitscn/Xsub0) ) 
END  IF 
floss  =  1.0 

c  If  target  type  is  'fluctuating'  calculate  fluctuation  loss  ratio. 

IF  (itgt  .EQ.  1) floss  =  ( (-LOG(psubd) ) * ( 1 . 0  +  gsubd/gsubfa) ) ** (-1) 
snmin  =  snmin  *  floss 

sndb  =  10.0  *  LOGIO (snmin) 

dbterm  =  10 . 0** ( ( 2 . O*antgn-rcvrnf-sndb-syslos) /lO . 0) 
rsubfs  =  58.0  *  (dbterm*pkpwr*sigma*tau/fmhz**2) ** . 25 
pthlos  =  32.45  +  20.0  *  LOGIC (rsubfs  *  fmhz) 
c 

c  Call  subroutine  prtout  to  print  out  input  values  and  calculated 

c  free-space  range  and  pathless, 

c 

call  rdrout (antgn, fmhz , hitsen, hbwidth, hsenrS , icoher, icmplx, itgt, 

1  pkpwr , plsrS , psubd , psubf a , rcvrnf , sigma , syslos , tau , 

2  rsubfs , pthlos) 

C  RETURN 

END 
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c  Subroutine  INRADR  prompts  the  user  for  radar  system  parameters 

c  and  returns  to  the  free-space  range  subroutine  FSRNG  to  deter- 

c  mine  the  free-space  range  of  the  radar  for  the  specified  in- 

c  puts  and  the  pathless  at  the  free-space  range, 

c 
c 

subroutine  inradr (antgn, fmhz,hitscn,hbwidth,hscnr8, icoher, icmplx, 
1  pkpwr,plsr8,rcvrnf ,syslos,tau) 

c 

real *4  antgn, fmhz , hit sen , hbwidth, hscnr8 , pkpwr , plsr8 , tau 
real*4  rcvrnf , syslos 
character*!  dummy 
integer*2  icoher , icmplx, ZW, ZR 
c 

c  Specify  the  read  (5)  and  write  (6)  channel  numbers. 

ZW  =  6 
ZR  =  5 

c  Enter  the  radar  system  parameters. 

write (ZW, '( "Enter  Radar  System  Parameters  ")') 
c  Set  the  radar  type  to  simple  radar  [  icmplx  =  0  ] . 

icmplx  =  0 
write (ZW, 1000) 

1000  format('Is  the  radar  a  simple  system  [range  and  azimuth  only,', 

1  /,' non-height-finder,  simple  signal  processing]  (y  or  n)?  ',$) 

read ( ZR , ' ( A1 ) ' ) dummy 
IF  (  dummy  .EQ.  'n'  )  icmplx  =  1 
c  Set  the  integrator  type  to  non-coherent  [  icoher  =  0  ] . 

icoher  =  0 

IF  (  icmplx  .EQ.  1  )  THEN 
write(ZW, 1005) 

1005  format ('Does  the  radar  use  non-coherent  integration  or', 

1  /, 'coherent  integration  [pulse  doppler  radar]  ' 

2  ' (n  or  c) ?  ' , $) 
read(ZR, ' (Al) ' )dummy 

IF  (dummy  .EQ.  'c')  icoher  =  1 
END  IF 
c 

write(ZW, 1010) 

1010  format (' Enter  frequency  in  MHz  (100  to  20,000)  ',$) 

read(ZR,*)  fmhz 

IF  (fmhz  .LT.  100.0)  fmhz  =  100.0 
IF  (fmhz  .GT.  20000.0)  fmhz  =  20000.0 
c 

write(ZW, 1015) 

1015  format (' Enter  peak  power  in  kw  (1.0  to  10,000)  ',$) 

read(ZR,*)  pkpwr 
IF  (pkpwr  .LT.  1.0)  pkpwr  =  1.0 
IF  (pkpwr  .GT.  10000.0)  pkpwr  =  10000.0 
c 

write (ZW, 1020) 

1020  f ormat (' Enter  pulse  width  in  microseconds  ', 

1  ' (>0.  to  1000.0)  ' ,$) 

read(ZR,*)  tau 
IF  (tau  .LE.  0.0)  tau  =  0.1 
IF  (tau  .GT.  1000.0)  tau  =  1000.0 
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c 

1025 

c 

1030 

c 

1035 

c 

1040 

C 

1045 

c 

1050 

c 

1055 


write(ZW, 1025) 

format (' Enter  receiver  noise  figure  in  dB  (>0.  to  20.0)  ',$) 

read(ZR,*)  rcvrnf 

IF  (rcvrnf  .LE.  0.0)  rcvrnf  =  0.1 
IF  (rcvrnf  .GT.  20.0)  rcvrnf  =  20.0 

write(ZW, 1030) 

format (' Enter  antenna  gain  in  dB  (>0.  to  50.0)  ',$) 

read(ZR,*)  antgn 
IF  (antgn  .LE.  0.0)  antgn  =  0.1 
IF  (antgn  .GT.  50.0)  antgn  =  50.0 

write(ZW, 1035) 

format (' Enter  assumed  system  losses  in  dB  (0.0  to  20.0)  ',$) 

read(ZR,*)  syslos 

IF  (syslos  .LE.  0.0)  syslos  =  0.0 
IF  (syslos  .GT.  20.0)  syslos  =  20.0 

IF  (icmplx  .EQ.  1)  THEN 
write (ZW, 1040) 

format (' Enter  the  number  of  hits  per  scan  (1  to  1000)  ',$) 

read(ZR,*)  hitscn 

IF  (hitscn  .LT.  1.0)  hitscn  =  1.0 
IF  (hitscn  .GT.  1000.0)  hitscn  =  1000.0 
ELSE 

write(ZW, 1045) 

format (' Enter  pulse  repetition  frequency  in  Hz 
' (1.0  to  100,000.0)  ' ,$) 

read(ZR,*)  plsrS 

IF  (plsrS  .LE.  0.0)  plsrS  =  0.1 
IF  (plsrS  .GT.  1.0E5)  plsrS  =  1.0E5 

write (ZW, 1050) 

format (' Enter  horizontal  beam  width  in  ', 

'degrees  (>0.0  to  90.0)  ',$) 

read(ZR,*)  hbwidth 

IF  (hbwidth  .LE.  0.0)  hbwidth  =  0.1 
IF  (hbwidth  .GT.  90.0)  hbwidth  =  90.0 

write (ZW, 1055) 

format (' Enter  horizontal  scan  rate  ', 

'in  rpm  (1.0  to  500.0)  ',$) 

read(ZR,*)  hscnrS 

IF  (hscnrS  .LT.  1.0)  hscnrS  =  1.0 
IF  (hscnrS  .GT.  500.0)  hscnrS  =  500.0 
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c 

c  Calculate  the  number  of  pulses  integrated  from  the  horiz.  beam 

c  width,  horiz.  scan  rate  and  the  pulse  repetition  rate, 

c 

hitscn  =  hbwidth*plsr8/6 . 0/hscnr8 
IF  (hitscn  .LT.  1.0)  hitscn  =  1.0 
END  IF 
c 

RETURN 

END 
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c 

c  Subroutine  INTRGT  prompts  the  user  for  target  parameters; 

c  target  size,  probability  of  detection,  probability  of  false 

c  alarms  and  if  the  target  is  fluctuating  or  steady.  It  then 

c  returns  to  the  free-space  range  subroutine  fsrng  to  deter- 

c  mine  the  free-space  range  of  the  radar  for  the  specified  in- 

c  puts  and  the  pathless  at  the  free-space  range, 

c 
c 

subroutine  intrgt ( itgt , psubd , psubf a , sigma ) 
c 

real*4  pexp, psubd, psubf a , sigma 
character*!  dummy 
integer*2  iexp, itgt , ZW, ZR 
c 

c  Specify  the  read  (5)  and  write  (6)  channel  numbers. 

ZW  =  6 

ZR  =  5 

c 

write(ZW, 1000) 

1000  format (' Enter  target  radar  cross  section  in  square  meters', 

1  '  (>0.0  to  l.e+07)  ' ,$) 

read(ZR,*)  sigma 

IF  (sigma  .LE.  0.0)  sigma  =  0.10 
IF  (sigma  .GT.  1.0E7)  sigma  =  1.0E7 
c 

write(ZW, 1005) 

1005  format (' Enter  probability  of  detection  (.1  to  .9)  ',$) 

read(ZR,*)  Psubd 
IF  (Psubd  .GT.  0.9)  Psubd  =0.9 
IF  (Psubd  .LT.  0.1)  Psubd  =0.1 
c 

write(ZW, 1010) 

1010  format (' Enter  probability  of  false  alarms  (4  to  12)  1.0e-',$) 
read(ZR,*)  pexp 
IF (pexp  .LT.  4.0) pexp  =  4.0 
IF(pexp  .GT.  12.0)pexp  =  12.0 
iexp  =  INT(pexp) 

Psubfa  =  10 . 0** (-iexp) 
c 

c  Set  the  target  type  to  'fluctuating'  [itgt  =  1]. 

itgt  =  1 
write(ZW, 1015) 

1015  format ('Is  the  target  fluctuating  [Swerling  Case  1]  or  ' 

1  /, 19X, ' steady  [Swerling  Case  0]  (1  or  0)?  ',$) 

read(ZR, ' (Al) ' ) dummy 
IF  (dummy  .EQ.  '0')  itgt  =  0 
c 

RETURN 

END 
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subroutine  rdrout (antgn, fmhz , hitscn , hbwidth, hscnrS , icoher , icmplx, 

1  itgt , pkpwr , plsrS , psubd , psubf a , rcvrnf , sigma , 

2  syslos,tau,rsubfs,pthlos) 
c 

c  Subroutine  RDROUT  provides  the  user  with  the  radar  system  para- 

c  meter  inputs,  the  target  size,  probability  of  detection  and  false 

c  alarm  rate  and  the  calculated  free-spaca  range  and  path-loss  values 

c  from  subroutine  RFSDR. 

c 
c 

real *4  antgn , fmhz , hitscn , hbwidth , hscnrS , pkpwr , plsrS 
real*4  psubd, psubfa , pthlos , rcvrnf, rsubfs , sigma , syslos , tau 
integer*2  icmplx, icoher , ZW, ZR 
c 

c  Specify  the  read  (5)  and  write  (6)  channel  numbers. 

ZW  =  6 
ZR  =  5 
c 

IF  (  icmplx  .EQ.  0)  THEN 
write(ZW, 1000) 

1000  format (///' Radar  is  a  simple  2D  system  (range  &  azimuth)  ') 

ELSE 

write (ZW, 1005) 

1005  format (///' Radar  is  not  a  simple  2D  system  ') 

END  IF 
c 

IF  (  icoher  .EQ.  1  )  THEN 
write(ZW, 1010) 

1010  format ( 'Radar  uses  coherent  integration  [pulse  doppler  radar]') 
END  IF 
c 

write (ZW, 1015)  fmhz 

1015  format (' Radar  frequency  =  ',f8.1,'  MHz') 
c 

write (ZW, 1020)  pkpwr 

1020  format (' Radar  peak  power  =  ',f8.1,'  kw') 
c 

write (ZW, 1025)  tau 

1025  format ( 'Radar  pulse  width  =  ',f6.1,'  microseconds  ') 
c 

write (ZW, 1030)  rcvrnf 

1030  format (' Radar  receiver  noise  figure  =  ',f4.1,'  dB') 
c 

write (ZW, 1035)  antgn 

1035  format ( 'Radar  antenna  gain  =  ',f4.1,'  dB') 
c 

write (ZW, 1040)  syslos 

1040  format ( 'Assumed  radar  system  losses  =  ',f4.1,'  dB') 
c 

IF  (icmplx  .EQ.  1)  THEN 
write (ZW, 1045)  hitscn 

1045  format (' Number  of  hits  per  scan  =  ',f6.1) 

ELSE 
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write (ZW, 1050)  plsrS 

1050  format (' Radar  pulse  repetition  frequency  =  ',f9.1,'  Hz') 
c 

write (ZW, 1055)  hbwidth 

1055  format (' Horizontal  beam  width  =  ',f4.1,'  degrees') 
c 

write (ZW, 1060)  hscnrS 

1060  format (' Horizontal  scan  rate  =  ',f5.1,'  rpm') 
c 

END  IF 
c 

IF  (  itgt  .EQ.  0  )  THEN 

write (ZW, '( "Target  is  steady,  non-fluctuating,  " 

1  "Swerling  Case  0")') 

ELSE 

write (ZW, '( "Target  is  fluctuating,  Swerling  Case  1")') 

END  IF 
c 

write (ZW, 1070)  sigma 

1070  format ( 'Target  radar  cross  section  =  ',pe7.1,'  square  meters') 
c 

write (ZW, 1075)  psubd 

1075  format (' Probability  of  detection  =  ',f3.1) 
c 

write (ZW, 1080)  psubfa 

1080  format (' Probability  of  false  alarms  =  ',pe8.1) 
c 

write (ZW, 1085)  rsubfs 

1085  format(//, 'Radar  free-space  range  =  ',fl2.1,'  km') 
c 

write (ZW, 1090)  pthlos 

1090  format (' One-way  path  loss  at  the  radar  free-space  range  =  ', 

1  fl2.1, '  dB' ,4(/)) 

RETURN 

END 
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Appendix  B 

FREE-SPACE  INTERCEPT  RANGE  (FSIR)  PROGRAM  LISTS 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


FSIR:  Calculates  free-space  ESM  intercept  range. 


VARIABLES : 


DESCRIPTION: 


antgn 

esmrng 

fmhz 

fsloss 

logtrm 

powrkw 

rsens 


Transmitting  antenna  gain  in  dB 
ESM  intercept  range 
Frequency  in  MHz 

Free-space  path  loss  at  esmrng,  in  dB 
Temporary  variable  -  log  of  powrkw  &  fmhz  terms 
Transmitter  output  power  in  kW 
ESM  receiver  sensitivity  in  dBm 


real*4  antgn,  esmrng,  fmhz,  fsloss,  logtrm,  powrkw,  rsens 


Call  inputs  to  enter  EM  system  parameters 


CALL  inputs (antgn , fmhz , powrkw, rsens) 


logtrm  =  10.0  *  ALOG10(  powrkw  /  (fmhz  *  fmhz)  ) 

esmrng  =  10.0**((  logtrm  +  antgn  -  rsens  +  27.5517  )  /  20.0  ) 

fsloss  =  20.0  *  ALOG10(  fmhz  *  esmrng  )  +  32.45 

Call  esmout  to  print  out  calculated  values  and  inputs 

CALL  esmout (antgn, esmrng, fmhz, fsloss, powrkw, rsens) 

stop 

end 
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subroutine  INPUTS (antgn,  fmhz,  powrkw,  rsens) 
c 

c  Subroutine  to  enter  the  necessary  inputs  to  determine  the  free- 

c  space  intercept  range  for  an  £SM  receiver, 

c 

c  VARIABLE  NAMES : 

c 

c  antgn 

c  fmhz 

c  powrkw 

c  rsens 

c  2R 

c  ZW 

c 
c 

REAL*4  antgn,  fmhz, 

INTEGER*2  ZW,2R 
c 

ZR  =  5 

ZW  =  6 

c 

WRITE  (ZW,1100) 

1100  FORMAT  ("This  program  will  calculate  ESM  free  space") 

WRITE  (ZW,1105) 

1105  FORMAT  ("intercept  range  and  path  loss  threshold.") 
c 

WRITE  (ZW,1200) 

1200  FORMAT  (/, 'Enter  frequency  in  MHz  (100  to  20000)  ',$) 

READ  (  ZR,  *  )  fmhz 

IF  (  fmhz  .LT.  100.0  )  fmhz  =  100.0 
IF  (  fmhz  .GT.  20000.0  )  fmhz  =  20000.0 
c 

WRITE  (ZW,1300) 

1300  FORMAT  (/, 'Enter  transmitter  power  in  kW  (0.1  to  10000)  ',$) 

READ  (  ZR,  *  )  powrkw 
IF  (  powrkw  .LT.  0.1)  powrkw  =  0.10 
IF  (  powrkw  .GT.  10000.0  )  powrkw  =  10000.0 
c 

WRITE  (ZW,1400) 

1400  FORMAT  (/, 'Enter  transmitter  antenna  gain  in  dB  (-5  to  50)  ',$) 

READ  (  ZR,  *  )  antgn 

IF  (  antgn  .LT.  -5.0  )  antgn  =  -5.0 
IF  (  antgn  .GT.  50.0  )  antgn  =  50.0 
c 

WRITE  (ZW,1500) 

1500  FORMAT  (/, 'Enter  receiver  sensitivity  in  dBm  (-10  to  -150)  ',$) 

READ  (  ZR,  *  )  rsens 

IF  (  rsens  .LT.  -150.0  )  rsens  =  -150.0 
IF  (  rsens  .GT.  -10.0  )  rsens  =  -10.0 
c 
c 

return 

end 


VARIABLE  DESCRIPTION: 

Transmitting  antenna  gain  in  dB. 
EM  system  frequency  in  MHz. 
Transmitted  power  in  KW. 

Receiver  sensitivity  in  dBm. 

Read  channel  number 
Write  channel  number 


powrkw,  rsens 


B-2 


subroutine  ESMOUT (antgn,  esiurng,  fmhz,  fsloss,  powrkw,  rsens) 
c 

c  Subroutine  to  enter  the  necessary  inputs  to  determine  the  free- 

c  space  intercept  range  for  an  ESM  receiver, 

c 

c  VARIABLE  NAMES:  VARIABLE  DESCRIPTION: 

c 

c  antgn  Transmitting  antenna  gain  in  dB. 

c  esmrng  ESM  intercept  range  in  km. 

c  fmhz  EM  system  frequency  in  MHz. 

c  fsloss  Free-space  path  loss  at  esmrng,  in  dB. 

c  powrkw  Transmitted  power  in  KW. 

c  rsens  Receiver  sensitivity  in  dBm. 

c  ZW  Write  channel  number 

c 
c 

REAL*4  antgn,  esmrng,  fmhz,  fsloss,  powrkw,  rsens 
INTEGER*2  ZW 
c 

c  Set  write  channel 

ZW  =  6 
c 

WRITE  (ZW,1000)  fmhz 

1000  FORMAT  (////,' Frequency  =  ',f8.1,'  MHz') 
c 

WRITE  (ZW,1100)  powrkw 

1100  FORMAT  ('Transmitter  power  =  ' , f 8 . 1 , '  kw') 
c 

WRITE  (ZW,1200)  antgn 

1200  FORMAT  ('Transmitter  antenna  gain  =  ',f5.1,'  dB  ') 
c 

WRITE  (ZW,1300)  rsens 

1300  FORMAT  ('Receiver  sensitivity  =  ',f7.1,'  dBm  ') 
c 

IF  (esmrng  .LE.  1.0e8)  THEN 
write  (ZW,1400)  esmrng 

1400  FORMAT  (/  , ' Free  space  range  =  ',F11.1,'  km') 

ELSE 

write  (ZW,1500)  esmrng 

1500  FORMAT  (/  , ' Free  space  range  =  ',lpel4.7,'  km') 

END  IF 
c 

write  (ZW,1600)  fsloss 

1600  FORMAT  ('Path  loss  threshold  =  ',F6.1,'  dB',5(/)) 
return 
end 
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Initial  Distribution 


Naval  Oceanographic  Office 
Stennia  Space  Center,  MS  S9622-5000 
Code  9200 


Naval  Enviroiunental  Prediction  Reaearch 
Facility 

Monterey,  CA  939i3-5006 


Center  for  Naval  Anatyaea 
Alexandria,  VA  22302-0268 


NOSC  Waahington  Liaiaon  Office 
Washington,  DC  20363-5100 


Defense  Technical  Information  Center 
Alexandria,  VA  22304-6145 


